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Abstract 

The electronic properties of many transition metal oxide systems require 
new ideas concerning the behaviour of electrons in solids for their explanation. 
A recent example, subsequent to that of cuprate superconductors, is of rare 
earth manganites doped with alkaline earths, namely Re\- x A x MnO^, which 
exhibit colossal magnetoresistance, metal insulator transition and many other 
poorly understood phenomena. Here we show that the strong Jahn Teller 
coupling between the twofold degenerate (d x 2_ y 2 and d Sz 2_ r 2) e g orbitals of 
Mn and lattice modes of vibration (of the oxygen octahedra surrounding the 
Mn ions) dynamically reorganizes the former into a set of states (which we label 
£) which are localized with large local lattice distortion and exponentially small 
intersite overlap, and another set (labelled b) which form a broad band. This 
hitherto unsuspected but microscopically inevitable coexistence of radically 
different I and b states, and their relative energies and occupation as influenced 
by doping x, temperature T, local Coulomb repulsion U etc., underlies the 
unique effects seen in manganites. We present results from strong correlation 
calculations using the dynamical mean-field theory which accord with a variety 
of observations in the orbital liquid regime (say, for 0.2 ~ x ~ 0.5). We outline 
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extensions to include intersite I coherence and spatial correlations/long range 
order. 

1 Introduction 

Solid state oxides containing transition metal ions with unfilled d shell electrons 
have been at the centre of attention in condensed matter and materials physics for 
the last two decades because of the variety and novelty of electronic phenomena in 
them and the possibility of new applications. The d electrons in solids are not as 
extended as s , p electrons or as localized as / electrons. Their motion in the system 
is highly constrained by the large local (on site) d — d repulsion or Mott-Hubbard 
correlation U. The central question is how this strong correlation, and other factors 
such as orbital degeneracy and electron lattice coupling, lead to electronic behaviour 
qualitatively different from that of conventional solids which are successfully described 
as degenerate Fermi liquids with well defined interacting electronic quasi-particles. 

While the best known examples of these are the high T c cuprate superconduc- 
tors [1 , during the last decade another group of oxides, namely the manganites 
Rei- x A x MnOz (where Re = La, Pr, Nd etc. and A = Ca, Sr, Ba etc.) has become a 
focus of major activity. The initial interest was sparked by the discovery [2] that their 
electrical resistivity changes enormously with the application of a magnetic field, the 
change (colossal magnetoresistance or CMR) being two or more orders of magnitude 
larger than the normal cyclotron orbital effect characterized by the dimensionless 
parameter (uj c t) 2 . Subsequent work has shown a bewildering variety of phases, phase 
transitions and phenomena [3J E| depending on the doping x, temperature T, and 
ionic species Re and A as well as external perturbations. An example is the phase 
diagram of Lai_ x C a x MnO% (Fig. [IJ which shows an insulating, Mn — O bond (Jahn- 
Teller) distorted but structurally ordered phase for small x, the transition of this at 
low T from a ferromagnetic insulator to a ferromagnetic metal at x ~ 0.2, and thence 
to a charge ordered insulating phase for x ~ 0.5. For 0.2 ~ x ~ 0.5 this becomes a 
paramagnetic insulator above a T c ~ 250 K, with CMR near T c . The phase diagram 
varies considerably with the ionic species. For example, La\^ x Sr x MnO^ has a param- 
agnetic metallic phase for 0.175 ~ x ~ 0.5 and shows no charge/orbital order , while 
Pr\- x Ca x Mn03 has no metallic phases (even ferromagnetic) [HJII|. Two other general 
characteristics are the following. First, physical properties are extremely sensitive to 
small perturbations; examples being the CMR itself, the unusually large strain and 
ion size effects 0011], the melting of charge/orbital ordering for anomalously small 
magnetic fields and a metal to insulator transition induced by the electronically be- 
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mgn substitution of 16 by O 18 0. Secondly, over a wide range of x and T, two very 
different types of regions, one insulating and locally lattice distorted and the other 
metallic and undistorted, coexist [5]. The regions can be static [TUl ITT] or dynamic 
[UJ[T21 US]; their size can vary from 100 A° [JJj to 3000 A° [10 . All these observations 
suggest that metallic and insulating phases are always very close in free energy 

These phenomena cannot be described satisfactorily in terms of well known mod- 
els for electrons in solids. For example, one might expect generally that doping 
LaMnO^ (a Mott insulator) with Ca introduces mobile holes in the Mn — d band 
so that the system ought to be metallic at least in the high symmetry paramagnetic 
phase; but it is not till x ~ 0.20. The observed properties of manganites need to be 
understood in terms of the active degrees of freedom which are believed to be the 
twofold degenerate e g electrons and the t 2g core-spins of Mn, and Jahn-Teller(JT) 
optical phonon modes of the oxygen octahedra and the interactions between them. 
There are three strong interactions present, namely the large on site d — d repulsion 
U ~ 5 eV [121 HE] amongst the e g electrons, strong JT mode - e g electron coupling 
g which splits the doubly degenerate e g level by about 2E JT ~ 1 eV [T7l[TH], and 
the large ferromagnetic Hund's rule coupling J H between e g and t 2g spins (~ 1 to 
2 eV jini El])- For comparison, the e g electron kinetic energy scale or band width 
W = 2D a is ~ 2eV |15| IT5| . The presence of a large JT coupling is clearly indi- 
cated by the two 'phase' coexistence [§] (with one 'phase' being lattice distorted and 
insulating, and the other metallic and undistorted) as well as by the sensitivity of 
properties to lattice strain or local disorder, and giant isotope effect jH] UH] alluded to 
above. Indeed, large JT distortions, while long ranged in LaMnO^, persist locally 

____ _____ rvj r^j 

well into the metallic regime [TT] IT2] (x < 0.3, and T > 77 K), where they are short 
ranged. Understanding the observed consequences of these interactions is one of the 
major challenges in the physics of strongly correlated electrons. 

Most existing theories [2H1 El 1221 E31 E] for manganites neglect one or the other of 
these strong interactions, make further approximations, and are inadequate qualita- 
tively and quantitatively. The earliest theoretical approaches [20J, commonly referred 
to as double exchange theories, consider solely the Hund's rule exchange Jh- How- 
ever, only a metallic state is possible in this case. A theory due to Millis, Mueller 
and Shraiman [21] additionally includes the effect of the coupling g, but treats the 
local JT lattice distortion classically, as annealed static disorder, and neglects U. 
A polaronic insulating phase also occurs in that case for large enough g, but the 
predicted results do not resemble experiments; for example at x ^ 0, one finds only 
metal-metal or insulator-insulator Curie transitions, unlike the commonly observed 
metal-insulator transition. The magnetoresistance is not colossal and there is no iso- 
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Figure 1: The phase diagram of La\- x C a x MnO^ in the doping x and temperature T 
plane (adapted from Ref. |24J). Various kinds of anti-ferromagnetic insulator (AF), 
paramagnetic insulator (PI), ferromagnetic insulator (FI), ferromagnetic metal(FM) 
and charge /orbit ally ordered insulator(CO/0 O*) are the phases shown. O and 
O* are Orthorombic(Jahn- Teller distorted) and Orthorombic(Octahedron rotated) 
structural phases. 
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tope effect. Dagotto and coworkers (see eg. Ref. j^j) have done extensive numerical 
simulations of several models, but on rather small lattices, and seen lots of instances of 
"phase separation". The early simulations explored the competition between (double 
exchange induced) ferromagnetism and anti-ferromagnetic (super)exchange, ampli- 
fication and generation of small scale phase co-existence by disorder, and identified 
the two kinds of magnetic domains with metallic and insulating regions. Later sim- 
ulations explored similar issues including static JT couplings. Based on these, they 
suggest jH] that transport in manganites should be pictured in terms of a random 
resistor network arising from tunnelling between misaligned ferromagnetic metallic 
domains across insulating regions, and that the CMR arises from enhanced tunnelling 
due to field induced alignment of the magnetic domains. There are many other mod- 
els designed to address specific effects, but no theoretical ideas which explain the 
novel general features of manganites cohesively. 

We propose here a new approach to the physics of manganites which includes 
all the three strong interactions, the double degeneracy of e g orbitals as well as 
the quantum dynamics of phonons. It leads to a simple physical description from 
which the properties of manganites follow naturally. Calculations of ground state 
and transport properties at finite temperatures are also presented. We believe that 
the ideas are of wide relevance to a large number of solid state and molecular systems 
with similar ingredients. 

2 Coexisting polaronic and band states 

Our ideas are principally based on two facts; namely the two fold degeneracy of the 
e g orbital and the strong (degeneracy lifting) JT interaction characterized by the 
large dimensionless ratio (Ejrp/hujo), (~ 10 for manganites, see ref.^7j) where hu>o 
is the JT optical phonon energy ~ 0.075 eV [25 . This coupling leads (as shown in 
greater detail below) to a large local JT distortion Q being associated with one linear 
combination of the orbitals ( JT polaron, labelled £ by us) when it occupied by a single 
electron. If the JT modes are not approximated as static displacements j2J but are 
treated quantum-dynamically, the inter-site hopping of the JT polaron is reduced 
by the exponential Huang Rhys factor 77 = exp{ — {E JT /2Tiujq)} ~ (1/200) for 
Ejt = 0.6 eV and huoo ~ 0.06 eV . This is the antiadiabatic limit overlap between 
the initial and final JT phonon wave functions (centred around Qo and 0) when 
the i electron moves in or out of a site. The corresponding effective ^-bandwidth, 
W* = ksT* ~ Wr] ~ fce(1251f), is thus very small, so the £ states are easily 
localized by any disorder present (eg. in cation site energy). They can hence be 
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regarded, to a first approximation (which may be inaccurate for T < T* because of 
inter-site coherence, see Section 8), as non dispersive localized levels. The £ polaron 
is probably close to this limit over a wide range of x and T, where large local JT 
distortion without long range order is seen ^TJ ^] . 

Since the e g orbital at each site is doubly degenerate initially, there is another 
orthogonal set of states which we label b, which have their largest amplitudes at 
the fraction x of hole sites where the polaron is not present; their occupancy on 
the polaron site costs a large extra energy U — (U + 2Ejt). The (bare) hopping i 
amongst these b states is not reduced and they form a broad band (of bare width 
2D ~ 2 eV ) whose properties are strongly affected, eg. their effective bandwidth 
2D is renormalized to smaller values, by the other two strong interactions present in 
the system, namely the repulsive scattering from the £ polarons (U) and the coupling 
to the t2 9 spins (Jh) depending on x and T. Roughly, 2D increases with x as well as 
with T _1 and H, because the inhibition of b hopping due to large U is reduced when 
there are more hole-sites, and that due to large Jh is reduced when the t2 9 spin order 
is enhanced. 

We believe that the unique feature of manganites is this necessary coexistence 
of antiadiabatic, JT distorted, localized (£) states and adiabatic, undistorted, broad 
band (6) states, arising from a spontaneous reorganization of the doubly degenerate e g 
states due to the large JT coupling and quantum phonon dynamics. This dynamically 
generated coexistence is qualitatively different from that of localized / and extended 
s,p, d electrons in rare earth solids, which has an atomic origin. We show below that 
for a wide range of x and T, both sets of states are partially occupied. This, we 
believe, underlies the ubiquitous two 'phase' coexistence observed [Hj. The high sen- 
sitivity of the physical properties of manganites to small perturbations arises because 
the (£) states, being localized and lattice distorted, are strongly influenced by local 
perturbations, which then affects the delicate relative stability of £ and b states. 

Finally, (again as discussed in more detail below,) the existence of the localized 
polaronic £ states in the presence of large U and Jh give rise to a new, major, 
doping dependent ferromagnetic nearest neighbour exchange coupling Jf between the 
t2g core spins. This comes about due to virtual, fast (adiabatic) hopping processes 
of the £ electrons to neighbouring empty sites and back, leading to a Jp roughly of 
order x(l — x)^ 2 /(2E JT S 2 ), the intermediate state energy due to the unrelaxed lattice 
distortion being 2Ejp. We get Jp — 2meV for x = 0.3. Our calculations suggest 
that this virtual, correlated double exchange due to the localized £ electrons, and not 
conventional double exchange due to mobile e g electrons as hitherto believed, is the 
dominant source of ferromagnetism and of the ferromagnetic transition temperature 
T c in the hole doped manganites. 
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3 A new model Hamiltonian for manganites in the 
strong electron lattice JT coupling regime 

Based on the above ideas, and for the purposes of making quantitative calculations 
and predictions, we have proposed [23 I2H] a new model Hamiltonian for manganites 
given by 

H a = £-(£jT + /i)C4,-/i£&, + A, (1) 

2,(7 2,cr 

- * £ b t b 3° + U £ n tiv n bia + H s 
{ij),a i,o- 

Here £f a creates the JT polaronic state of energy —E JT and spin cr, localized at site 
i (£ ia is the corresponding destruction operator). The broad band electron, created 
by the operator bf a , has mean energy zero, and nearest neighbour effective hopping 
amplitude t. The two repel each other on site with energy U. The common chemical 
potential /i is determined by the constraint that the filling, i.e., the average number 
of e g electrons per site must be determined by the doping x according to 

N- 1 -£(m) = N- 1 £ £((n to ) + (nua)) = (1 - x). (2) 

i i O 

The term H s in Eq. El models the spin dependent interactions, and is given by 
H s = -Jn^Si ■ Si - J F £ Si ■ Sj - [iB^Si ■ h (3) 

i <ij> 

It includes the strong ferromagnetic Hund's rule coupling Jh between the e g spins 
Si(= (su + Sbi)) and the (spin S = |) t 2g spins Si , the net effective ferromagnetic 
exchange coupling J F between these t 2g spins (from the mechanism alluded to above, 
and discussed in more detail below), and the interaction of the latter with an external 
magnetic field H |2H]- For simplicity, in what follows we further approximate the t 2g 
spins as classical fixed length spins whose directions fluctuate, i.e. write Si = Sfli 
where f2j is a unit vector. 

We now discuss how the above model Hamiltonian can be motivated from the 
microscopic Hamiltonian for the manganites in the limit of strong JT interaction. 
A reasonably realistic microscopic model Hamiltonian or energy operator H of d 
electrons (e g and t 2g ) in the manganites has three types of contributions, namely 
those involving only the e g electrons (H e ), the coupling of these electrons with the 
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JT lattice modes as well as the energy of these lattice modes themselves (H e i), and 
the part H s which involves the e g and t 2g spins. Thus 



H = H e + H el + H s (4) 

where 

H e = J2( 6 i ~fj)nia+ 3 £*ii S £ ( 5 ) 

i,cr <ij> 

+ U (fkafkv' - ^io) 
i,<T<rf 

H a = '£(H i JT + Hi) = gJ2(at T T z a ia Q M + at a T x a ia Q 2 ^ (6) 

i i,cr 

i i 

Here the operators af a , di a respectively add and remove an electron at site i, in the 
spin state a which can take two values, and in the two e g orbital states with labels 
a = 1 corresponding to d 3z 2_ r 2 and a = 2 corresponding to d x 2_ y 2] i.e., af a is a short 
hand for (af la , af 2a )- The number operator n ia is the usual sum ^ Q af aa a iaa = af a a ia . 
The e g electron has site energy e$, equal to for a clean system. The (anisotropic 
and orbital dependent) nearest neighbour hopping matrix is and U is the repulsion 
between electrons at i in different states. The first set of terms in H e £ describes the 
coupling of e g electrons with the two JT lattice modes Q 3i and Q 2 i with strength g, 
[r a being the usual Pauli matrices in orbital space) the former displacement leading 
to a splitting and the latter to a mixing of the two e g states. The last two terms in 
H e i are the potential and kinetic energies of the modes neglecting anharmonic and 
intersite terms. H s is the same magnetic Hamiltonian as in eqn. El (except that, 
in the context of the microscopic hamiltonian, the e g spin Sj = \ X^// ^(^W'^m' 
where a is the Pauli spin operator). 

In manganites the JT coupling g is large compared to t/Q Q (where Qq is the typical 
size of the JT distortion), and one can first concentrate on the single site, coupled 
electron-phonon problem. The Schrodinger equation corresponding to H l JT + H\ can 
be solved exactly at a single site. Qi = (Q 2 i, Qsi) is like a local pseudo- magnetic field 
splitting the pseudo-spin 1/2 orbital levels 1 and 2, hence H l JT has eigenvalues ±gQi, 
where Qi = \Qi\ is the magnitude of the JT distortion. The diagonalization of H l JT 
can be done explicitly using the operators 

b+ = cos{e i /2)a ¥ il + stn(6i/2)a + l2 (7) 
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and 

If = -sin{Bi/2)aUx + cos{6i/2)a + i2 , (8) 

where 9i = tan~ 1 (Q 2 i/Q3i), i.e., it determines the orientation of the JT distortion. 
When an electron is present at site i in the £ state, the effective lattice potential 
energy at site i in H JT + H\ has the form 

Vr = (KQ 2 J2 - gQi) (9) 

This has a minimum at Qq = (g/K) (the JT polaronic distortion) with an energy 
lowering of Ejt = (g 2 /2K) ~ .5 eV independent of Qi jSUj. The JT polaronic £ state 
that we referred to above is this £ state together with the JT distortion. When an 
electron is present in the other 'anti' JT , b state, the effective lattice potential energy 
is 

V l + = (KQ 2 /2 + gQ l ) (10) 

with a minimum at zero displacement so that no distortion is associated with it. 

Thus, when an I polaron is present at a site, low energy lattice configurations 
around the deep minimum of V[~ involve a shift of Qi by an amount Qo at each site. 
This is formally achieved by a Lang-Firsov transformation 31\ fji = exp { — {ipiQo/h)} 
where Pi is the radial momentum operator conjugate to Qi. In other words, formally, 
one has if = fjilf . The Huang-Rhys reduction factor |26] r] = exp { — {Ejt /2TujJq)} 
we have alluded to earlier, which multiplies the hopping involving the £ states, is just 
the ground state expectation value of fji, corresponding to the anti-adiabatic limit 
overlap of the distorted phonon wavefunction at a site with the undistorted phonon 
wavefunction. (We ignore fluctuations on the basis that irj -C huj [32] )• The coherent 
intersite hopping of the £ electrons is thus t* — if] ~ 10 K, very small, and can be 
neglected to a first approximation. The 'anti' Jahn Teller state b has no associated 
lattice distortion, so its hopping amplitude Uj is not reduced, but depends on the 
angles ^ and 9j at sites % and j. We assume a nearest neighbour hopping amplitude 
iij = i averaged over angles ^ either statistically (and classically) or by quantum 
fluctuations . Such an orbital liquid approximation is reasonable for the phases with 
no long range orbital order, i.e., for 0.2 ~ x ~ 0.5 in most manganites, but poor for 
those other values of x for which one has strong orbital correlations or long range 
orbital order. 

The localized polaronic £ and the broad band b states coexist in manganites 
since the conditions for anti-adiabaticity of the former, namely it] Hlvq and the 
adiabaticity of the latter, namely i 3> huo, are both fulfilled. The b electrons are 
repelled strongly by the £ electrons, with energy U = U + 2E JT , and hence will have 
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the largest amplitudes at sites where the £ electrons are not present, i.e. at the hole 
sites. They hop quickly (time scale h/t <C l/u) Q ) among such sites, avoiding strongly 
the sites with static £ electrons, where they hence have small amplitudes [33 . 

Finally we discuss the result alluded to above namely that the existence of local- 
ized, JT distorted £ states in the presence of large Jh and large U can give rise to a 
new mechanism of ferromagnetic exchange between the t 2g core spins, which we refer 
to as "virtual double exchange". 

Suppose that a localized i electron is present at a site, and that site is distorted. 
The £ electron can take part in fast (adiabatic) virtual hopping processes to neigh- 
bouring sites, i.e., leaving the local lattice distortion unrelaxed, by paying an energy 
cost of 2Ejt in the intermediate state. For large U this can happen only if the 
neighboring state is empty, and for large J H only if the t 2g spins on the two sites are 
parallel. (Otherwise the energy of the intermediate state will increase by U and Jh 
respectively.) Clearly, from second order perturbation theory this process will give 
rise to a new, £ occupancy dependent, ferromagnetic exchange coupling between the 
t 2g core spins of the form 

( P \ 1 * - 

-(Oj • Qj + l)[n a (l - rij) + n ej (l - n*)] (11) 



\2E JT J 2 

(where we have again ignored the dependence on the angles 9t etc). The |(fij ■ Clj + 
1) factor comes from large Jh, and the occupancy dependent terms from large U. 
Within the homogeneous orbital liquid approximation we are using in this paper, this 
translates to an effective ferromagnetic interaction proportional to xfii ~ x(l — x) as 
stated earlier. 

The normal super-exchange coupling between t 2g core spins on neighboring sites is 
of order P/U , and its sign depends on the nearest neighbour orbital correlations [3*lj . 
For example, in LaMnO^ it is anti-ferromagnetic along the c axis, and ferromagnetic 
in plane. Conventional double exchange can, in principle, explain metallic ferromag- 
netism in the manganites, but would be unable to account for ferromagnetism in 
their insulating phases in the absence of orbital correlations. The above mechanism 
yields a ferromagnetic coupling which is clearly much larger than and dominates the 
super-exchange for all intermediate doping, operates even in the insulating states 
and even in the orbital liquid phase, and as we discuss later, makes the dominant 
contribution to the ferromagnetic T c even in the metallic case. 

Putting all this together, expressing the on site coulomb and exchange interactions 
in the £ and b basis, and including the new exchange mechanism, leads us to the 
suggestion that in the large {Ejt/TioJo) limit, and for describing low energy electronic 



10 



properties the microscopic model Hamiltonian in Eq. 0] can be approximated by the 
much simpler, two fermion species, Hamiltonian Ha, given in Eq. El The relatively 
small gf^BH ■ Sj term and other coulomb repulsion terms (between I and b electrons 
of different spins) have been neglected in Ha, because of the large Jh (Jh ^>t) limit 
of relevance here (whence the spin of the I and b electrons is forced to be parallel to 
the t 2g spin), and we have relabelled the effective repulsion between the two types of 
electrons as U. 

4 DMFT wth polaronic and band states 

Apart from the constraint equation Eq. |2] and the spin dependent H s , Ha, is es- 
sentially the well known Falicov-Kimball model (FKM)jnS] for non hybridizing / 
electrons in a broad band metal with correlations. Indeed at T = 0, in the ferro- 
magnetic phase and for large Jh, the spin degrees of freedom are completely frozen, 
whence it reduces to the (spin-less) FKM. More generally, it describes the dynamics 
of the b electrons moving in an annealed disordered background of immobile I elec- 
trons and t 2g spins in the presence of strong on site repulsion U and Hund's coupling 
Jh , the annealed disorder distribution being thermodynamically and self- consistently 
determined. This strong correlation problem cannot be solved in general, but can be 
solved within the framework of dynamical mean field theory (DMFT) |36j which is 
exact at d — oo, and is quite accurate for three dimensions. We have carried out 
these DMFT calculations for thermodynamic properties (the occupancies fie, fib, the 
magnetization m = (Si), the specific heat, magnetic susceptibility, etc.), spectral 
behaviour (eg. the b band self-energy, propagator and density of states(DOS)) and 
transport (resistivity and magnetoresistance). While the specific results we discuss 
below come from such calculations, we believe that many of the physical mechanisms 
that underlie the results have a larger sphere of validity. 

The properties of Ha, are determined largely by the dynamics of the b electrons, 
eg., their propagator or Green's function Gij(ui), as affected by the (annealed dis- 
ordered) distribution of n& and Oj due to strong U and Jh- In the DMFT frame- 
work, the b electron self energy due to these interactions is site local, i.e. 
= <5ijXXcu), and is determined self consistently. The general approach is well 

known 

Essentially, the problem reduces to that of one site, with its local degrees of free- 
dom and interactions, immersed in an effective medium (comprising the other sites). 
We make the simplest reasonable approximation corresponding to a homogeneous 
annealed system, namely that (n&j) = , (nua) = ^&<r and (Oj) = m = mz are the 
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same at every site i, and that there are no correlations between these quantities at 
different sites jH7]. Consistent with this, we approximate the ferromagnetic exchange 
interactions in H s in terms of a homogeneous molecular field in the standard way 
and replace 

H s = -J H SY,Si J F m) -J2&i (12) 

i i 

Here h = gfisSH , Jp = IzJpS 2 where z is the co-ordination number. The quantity 
m, proportional to the magnetization, and ng a , the average occupancy of the polaronic 
states, are also to be determined self-consistently. 

Thus, in the present context, the effective medium is a homogeneous electron bath 
with which the local b electron hybridizes, the molecular field Jpm in Eq. El and the 
local distributions of ri£ a and fl . The on-site b electron propagator due to the medium, 
but without the local interactions, is Q where the (Matsubara) frequency variable 
and the spin label are suppressed. The single site problem can be solved exactly, just 
as in the Falicov-Kimball case [23, to determine (1) the (homogeneous) annealed 
distributions Pint) and -P(O) as functionals of Q, m and the model parameters and 
(2) the local b electron propagator G as a functional of Q ,m, the model parameters 
and these annealed distributions. There are two other formal relations among G, 
G and E, namely (G" 1 = Q~ x - S) and G = £g[(G£ -1 - By iterating these 

equations to self-consistency , one can determine S, G, Q , P{nn) and P(Q) (and hence 
nta, and m) explicitly for any chosen values of the model parameters U, eg, t and 
Jf at a fixed T and \x (the latter being chosen such that ne + flb = (1 — x)). Thermal, 
spectral and transport properties can be calculated in a straightforward way [UHl from 
the b electron propagator. For example, the current-current correlation function (the 
Kubo formula) which determines the electrical (and optical) conductivity, can be 
expressed entirely in terms of Gij(u), vertex corrections being negligible. 

The resulting calculations, while straightforward, still involve extensive numerics. 
The calculations whose results are described below have been done in the large Jh 
limit, i.e. assuming that the spin of the I and b electrons are forced to be aligned 
along t2 9 spin Q, whence there is considerable simplification [?Ill27j. We have done 
such calculations both for a realistic model t and for a model semicircular density of 
states. The two results are not very different, and here we mostly discuss results for 
the latter case, which has the advantage that some exact analytical results can be 
obtained in the U = oo, T = limit, and even the finite U, non-zero T calculations 
are much simpler. The details are described in refs. |2H1 12*7| 135] and other papers we 
will publish elsewhere. 
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5 Metal insulator transitions 



To start with, we show how the ideas above lead to a simple physical understanding 
of insulator metal transitions (IMT) in manganites, based on the variation of the 
occupancies and relative energies of the I and b states with U, x, T etc. 

Fig. El shows the density of states for (Ejt/D ) = 0.5 , D Q = 1 eV , U — oo (and 
Jp = 60.2 meV) for several values of x at T = 0. At T = 0, because of Jf the ground 
state is ferromagnetic, and the spin degrees of freedom are completely frozen. For 
large U, the renormalized or effective b bandwidth D self consistently goes to zero 
as x — > 0, whence for small x the bottom of the b band will be above the I level 
(of energy -Ejt) and all electrons are in the latter states, as in Fig. 2a. For, the 
low energy b band states are constrained to reside mostly on the small fraction x of 
empty, undistorted sites, and have a small amplitude at the lattice distorted sites 
(where I electrons are present) from which they are repelled with energy U . The 
system is hence an insulator. This cannot happen in a pure double exchange model, 
where ferromagnetism and metallicity go together. As x increases, so does D, and 
there is a critical x c (= 0.25 within DMFT for the parameters used in Fig. 2) at which 
the b band bottom touches the I level energy as in Fig. 2b. For x > x c , the b band 
bottom goes below the I level, and some of the I electrons are transferred to the b 
band until all b levels up to the (now renormalized) I level are occupied as in Fig. 
2c. The system is then a ferromagnetic metal, but most of the electrons (0.6 per site 
out of the 0.7 present) are in the polaronic i state though there is no long range JT 
order. 

Thus there is an IMT, not at x = as in a naive doped Mott insulator picture , 
but at x c 0. This shift is due to the strong electron lattice coupling g which makes 
the JT polaronic £ level possible and stabilizes it by an energy E JT , and large U, 
which makes the b electron bandwidth vanish as x — > 0. For U — > oo, within DMFT 
we can show [28] that D = D ^fx whence we obtain x c analytically as x c = (Ej T /D ) 2 
where D Q is the bare half-width of the b-band. The decrease of x c with increasing 
D Q and its numerical value for reasonable choices of Ejt and D Q are in accord with 
experiments (for example, as one goes from LaCaMnO^ to LaSrMnO^ which has a 
larger bandwidth). Finite U corrections to x c , which we have calculated, are small 
for U ~ 5eV which is appropriate to these systems. A more detailed discussion of 
these and other zero temperature results is presented in [2~8] . 

FigEl shows the evolution of the spectral functions with temperature, for parame- 
ter values that correspond to Sr doping, for a fixed doping of x = 0.175 > x c for this 
system, and for the temperature values T =100K, 180K, 208K, 300K. For x > x c , the 
system is metallic so that both b and I electronic states occupied. On increasing the 
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Figure 2: Evolution of DOS of the b band for various doping values of x — 0.1, 0.25, 
0.3 and T = with parameters E JT = 0.5 eV, D D = leV. The effective £ level, 
labelled as e% is also indicated. 
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spond to Sr doping (see text), for T = 100K, 180K, 208K, 300K and x = 0.175. 
Thick lines represent the effective £ level and up and down arrow indicate up spin 
and down spin spectral functions. 
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temperature, the t 2g spin order decreases and the b band narrows again, since due 
to double exchange the effective hopping amplitude of the b electrons decreases. 
The spectral weight for the down spin polarization also turns on, reflecting the reduc- 
tion in the magnetization from its zero temperature value. The sequence of figures 
3a to 3c illustrates this progression. Eventually, both up and down spin spectral 
functions become equal, and the b band bottom crosses the I level and moves up as T 
rises beyond T c (= 208 K in the present context, see Fig. 3c) where ferromagnetism 
disappears, and the system becomes a paramagnetic insulator as in Fig. 3d, with b 
states occupied thermally across a gap. We thus have a simple picture of the thermal 
ferromagnetic-metal to paramagnetic-insulator transition. 

6 Resistivity, CMR and material systematics 

There are very few calculations El 122] of the transport properties of manganites 
though this is one of their main unusual features. Within the DMFT for our model, 
the electrical resistivity p can be calculated in a straightforward wav|36| from the 
b electron propagator or Green's function. Fig. 4 shows some of our results, for 
x = 0.3 and parameters representative of Nd\- x Sr x MnO^ and La\- x Ca x MnO^ [4*Uj . 
Experimental values are also given for comparison |41j . 

We see that the resistivity of the paramagnetic insulating state is fairly well de- 
scribed by the theory. In this phase, the conducting b band is only thermally occupied. 
The effective electrical gap ~ 850-fT (the experimental value is ~ 12507^01]). The 
positive feedback between the thermally excited, strongly T dependent number n^iT) 
of b electrons, and the double exchange broadening of the b band as T decreases, leads 
to a very rapid decrease and closure of the electrical gap just below T c and the con- 
sequent rapid decrease in the resistivity just below T c as in experiment to a rather 
large value ~ 2 mVtcm. But we note that the calculated resistivity does not decrease 
much thereafter down to T = 0, unlike experiments, where p{T) decreases to residual 
values of ~ 50 pQcm below about T = 125K. We believe that the latter is due to 
inter-site i coherence, neglected here (see later). 

We also show in the Fig. 0]the resistivity in a field of ITesla for Ndi_ x Sr x MnO^. 
We clearly see CMR. The external magnetic field polarizes the t 2g spins, and via Ju 
increases the b band width which reduces the electrical gap. CMR arises because 
resistivity depends exponentially on the gap. The effect is largest near T c as the 
change in the magnetization due to a given magnetic field is the largest there. As our 
mean field approximation for the Curie transition neglects short range magnetic order 
above T c , p is somewhat overestimated and the CMR somewhat underestimated. 
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Figure 4: Electrical resistivity and CMR in the coexisting JT polaron - broad band 
model for Ejj- = 0.5 eV, U = 5eV and x = 0.3. D a and Jp are chosen so as to 
reproduce the experimental |H] T c and p(T c ) of Ndi_ x Sr x Mn0 3 (D a = 1.05 eV, Jp = 
52.7 meV) and of La^Ca^MnOs {D = 1.15 eV, J F = 60.2 meV) . Calculated p(T) 
at H = ITesla is shown for Ndo^Sro^MnOs 
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The properties of manganites vary strongly and characteristically with the ionic 
species Re and A. The material systematics of such variation, eg., in the thermal 
IMT as well as the CMR depend in our model largely and sensitively on the ratio 
(D /Ejt) and somewhat on the D a and Jf- A small increase in D q /Ejt reduces the 
high temperature resistivity p(T > T c ) enormously, as is clear in Fig. 0] where we 
have increased D Q from 1.05 to 1.15 eV and Jp from 52.7 to 60.2 meV to fit T c and 
p{T c ) for La 7 Ca , 3 MnO 3 . Broadly this is because the density of current carrying 
b electrons and hence cr(T) = l/p(T) depends exponentially on (Ejt — D Q ). This 
is also the general reason for the observed large variation of p(T > T c ) and T c with 
strain and pressure. 

Making the reasonable assumption that Ejt, a single octahedron quantity, is 
roughly unchanged across materials, but that Jp, the effective ferromagnetic ex- 
change, scales as D%, we have calculated properties as a function of D Q for x ~ 0.3. 
We find that T c is typically fairly close to the (Curie- Weiss) value (Jf/3) (and hence 
increases as D 2 . as D Q increases), with a slight enhancement due to a nonzero nf,(T c ) 
and double exchange. With increasing D , the resistive transition on lowering T 
through T c changes from insulator-insulator to insulator-metal, and then to metal- 
metal. These trends are seen experimentally [6, as a function of increasing cation 
radius (ta), which is known to roughly track the bandwidth D Q . We also find that 
the relative CMR, defined as [R(H) — R]/R(H), depends exponentially on T c (in- 
creasing as T c decreases). This is indeed the behaviour extracted from measurements 
on several systems |%2] . The observed unusually large dependence of T c and of the 
resistivity at T c on pressure (12] in Lai_ x Ca x Mn0 3 is reproduced if we assume that 
D increases with pressure at a realistic rate of about .01 eV/Kbar. A more detailed 
discussion of these finite temperature results and material trends is presented in |27j . 

A key feature of our theory is the result mentioned above that the dominant 
contribution to T c comes from virtual double exchange of the localized I polarons 
and not from conventional double exchange, since the mobile carrier density (n^) 
is very small. Interestingly, the precipitous, universal and nearly linear drop in T c 
with local ion size variance reported by Attfield[7] in a variety of multi-component 
manganites is much larger than what is expected from a simple double exchange 
picture. In the double exchange mechanism, T c is proportional to which depends 
on the Mn — O — Mn bond angle 0^-, which is strongly affected by the local ion-sizes, 
as cos(4>ij). The variance in 0^-, namely < (5<fiij) 2 > x l 2 is at most 10 degrees, so 
that one expects that the resulting reduction in T c is AT C /T C ~< (8(fiij) 2 > /2 ~ 
1/60. The observed reduction is an order of magnitude larger! We believe that 
such a large reduction is however understandable within our theory if we extend the 
calculations described above to include the effect of ion size variance on Ejt , and 
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£ — b hybridization effects as discussed in section 8. 

7 Other unusual properties 

We now illustrate through three examples how other unusual properties of manganites 
can be understood in our approach. Specifically, we consider the anomalously low 
carrier density as inferred from optical conductivity in the metallic phase [44J, the 
small electronic specific heat jl] and the electron hole asymmetry [4*3] . 

7.1 Low metallic carrier density 

In the metallic regime (0.2 ~ x ~ 0.5), one can estimate n e ff , the effective density of 
carriers, from the the Drude weight i.e. the frequency integral of the real part of the 
optical conductivity a r (uj). In a free electron like limit, this integral is (rcne 2 /2m*) 
where n is the density of carriers, and m* is the band optical mass. Since there is 
no evidence for enhancement of m* with respect to the band mass (eg. see Ref. 
results for La 1 ^ r Sr T MnO ? . [44] lead to n e ff = 0.06 for x = 0.3 and T = 0. The 
expected value (for a doped Mott insulator) is (1 — x) = 0.7. Although there are 
several uncertainties in inferring n e // precisely from a r (uj), its reduction by nearly an 
order of magnitude, and even more so its observed strong decrease with temperature 
on the scale of T c which is much less than the (normal) Fermi temperature, are 
puzzling. But the results are easily understood in our model, since the low frequency 
cx(u;) is due to b electrons. As seen from Fig. (2c) at T = there are indeed 
very few electrons in the b band, all near its bottom, below the £ (or Fermi) level. 
The DMFT calculation with parameters appropriate for Lai- x Sr x MnO^ at x = 0.3, 
namely E JT = 0.45 eV , D Q = 1.28 eV , J F = 1150 K (and consequently T c = 390 K), 
leads to n h {T = 0) = 0.12, and a decrease with temperature which is very similar to 
what is seen in experiment. 

7.2 Normal linear electronic specific heat 

The low T electronic specific heat in manganites (see Ref. |4] for a review) goes as 
7T with a 7 nearly independent of x, unlike other doped Mott insulators(eg Ref. [I]) 
and has a value expected for a tight binding band with bandwidth 2D Q ~ 2.0eV. In 
polaron theories with fii polarons per site C v = Ughs in the classical gas limit, and 
~ nik B (k B T/D*), where D* is the effective bandwidth, for T (D*/k B ). Neither 
of these is observed and this is often argued to be evidence against small polarons in 
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manganites. In our model, although the electrons are mostly JT polarons, we expect 
that they will make very little contribution to the low temperature specific heat. 
For, being localized they will strongly couple to disorder potentials and affected by 
coulomb interactions, and we expect that their entropy will be frozen out at a large 
temperature scale. We find that 7 (calculated in the paramagnetic metallic phase to 
avoid large magnetic contributions to the specific heat) is about 0.9 (m J/ K 2 mole) 
with the parameters used above for La 7 Sr 3 MnO 3 . The value for a band of width 
2eV is about 2.0 (mJ/K 2 mole). The physical reason for the smaller 7 is that the 
few b electrons which are present occupy states near the bottom of the b band with 
energy close to —Ejt, and the density of states there is relatively smaller (Fig. 2c). 

7.3 Electron hole asymmetry 

The electronic and magnetic properties of hole doped (2 < 0.5) and "electron doped" 
((x > 0.5) manganites are very different. For example, the latter are much more 
metallic than the former in the paramagnetic phase [15] (although the AF or charge 
ordered phases can be insulating for other reasons). This behaviour is unexpected 
in a one orbital strong coupling antiadiabatic polaron model since for low e g carrier 
density the dilute small polaron assembly will form an insulator, as also in a model 
with only Jh , where there is electron hole symmetry. It is often stated that polaronic 
effects are not seen for e g orbitals in the e doped regime because the JT coupling 
becomes weaker as x — > 1. This seems unlikely given the local nature of the JT 
interaction and the stability of the (MhOq) octahedron. 

The metallic nature of the e doped regime which has low e g density is natural in 
our theory because for large x the effectively uncorrelated b electrons form a wide 
band whose bottom is occupied by the small number of e g electrons, and the I states 
are unoccupied, as depicted in Fig. 2d. This fact also explains why pure band models 
are fairly successful [46] in describing the magnetic ground states in the the e doped 
limit. 

8 Discussion 

In summary, we have presented here a new model of coexisting localized JT polarons 
and broad band electrons for manganites, argued that it arises inevitably in the 
presence of orbital degeneracy and strong JT coupling, and shown that it explains a 
wide variety of characteristic properties of manganites. In this concluding section we 
compare our theory with some of the other theories, discuss some of the inadequacies 
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of the present model and extensions to overcome these. We believe that a more 
general treatment of the model Eq. |2]with some extensions can lead to a complete 
description of manganites. Two examples, namely inclusion of spatial correlations 
and of intersite £ state coherence are discussed below. 

To begin with, we note that the picture developed here is very different from 
that in other polaronic models, which either neglect e g orbital degeneracy [221 12S] or 
work in the adiabatic approximation |2*T) 147] . In the former, for example, at high 
temperatures transport is due to the activated hopping of localized small polarons. 
In the latter (adiabatic) models, the polarons also form a broad band, whence it is 
difficult to obtain a paramagnetic insulator to ferromagnetic metal transition say at 
x ~ 0.3. Furthermore, small polarons disappear below T c even for large g, the Drude 
weight is not small, and there is no isotope effect, all in disagreement with experiment. 
In both, it is argued that small polarons are likely only at high temperatures, and 
that small and large polarons coexist only in a narrow range of x and T determined 
by the effective electron phonon coupling and that there are no small polarons at low 
T. 

In our theory, the carriers are broad band electrons thermally promoted out of 
localized JT polaronic states. Small JT polarons (in the anti-adiabatic limit, with 
negligible bandwidth) and band states(in the adiabatic limit) necessarily coexist over 
a wide range of x (0.2 < x < 0.7 say) and for all T. In contrast, in single orbital 
polaron models, there is a crossover with increasing dimensionless electron phonon 
coupling X e ff(x,T) = {g 2 / K W e ff(x,T)} from large to small polarons, so that the 
two can occur together only over a narrow range of x and T. 

An unusual feature of the polaron level in our model is that it has no prominent 
thermodynamic or spectroscopic (sharp level) manifestations. Firstly as pointed out 
above, it does not lead to a large electronic specific heat as we expect that their 
entropy will be frozen out at energy scales determined by coulomb and impurity 
interactions not explicitly included in our model. The remaining linear specific heat 
for the b band metal, is roughly of the size seen experimentally^. Secondly, no 
sharp £ like spectral feature is to be expected. Rather, the £ excitation spectrum 
will be an incoherent continuum, starting from a weak (weight oc rj ~ 1/200) anti- 
adiabatic low energy part building up to adiabatic (Franck-Condon) higher energy 
features, as the fast removal of an £ electron leads to highly excited lattice states 
(energy ~ 2E JT ). Indeed, optical conductivity data in manganites do showjH] such 
an incoherent continuum which is largely independent of temperature, and we believe 
that this is their origin. 

In our theory presented earlier we made the simple approximation that the local 
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JT distortion Q Q is independent of x and T, and that its large value is determined 
entirely by on-site properties, namely the electron JT phonon coupling g and the 
phonon mode force constant K. This is probably valid for high temperatures not too 
far below T c and for x < 0.5. However, JT polaronic effects can become dynamical 
and unobservable at low T for two reasons . Firstly, the local lattice kinetic energy 
has a term (1/2)MQ 2 9 2 where Q D is the magnitude, and 9{ the direction or orienta- 
tion, of the local JT distortion. Hence which also determines the orbital mixing 
amplitudes of the t and b states, rotates due to quantum fluctuations, although due 
to anharmonicity and crystallinity, the rotation could be hindered and slow. A sec- 
ond, perhaps more important, source of this dynamical J — T effect is the quantum 
mechanical intersite I state coherence which makes the 9i ill defined. The energy 
scale for the latter, D* ~ 125K, is consistent with the widely observed |H I§1 ITT1 IT3] 
reduction or disappearance of static or of long time scale local JT distortion in metal- 
lic manganites as their temperature decreases well below 125K. Finally, intersite i 
coherence, especially hybridization with energetically degenerate broad band b states 
in the metallic regime, will reduce Q Q and broaden the I band. The two effects feed 
back on each other, and it is in principle possible that small polarons may weaken and 
disappear altogether at low temperatures in the metallic phase rather than merely 
becoming dynamic. We have not explored this possibility here. 

Intersite i state coherence can be included in our model by adding to Eq.(l) an 
t — b hybridization term H(y),<r ti(£ + ia bj a + b + j a £ ia ) + t 2 £f a £j a where t\ ~ irj ~ D*, 
and t 2 ~ irj 2 . We expect that its inclusion, and the consequent development of long 
range intersite I state coherence at low temperatures will, in addition, lead to the 
observed smooth decrease [19J in the electrical resistivity of clean metallic manganites 
from about 1 mflcm just below T c to a small value, ~ 50 — 100 fiQcm, as T — > past 
a characteristic crossover temperature of about 100K, of the same scale as D*. In the 
DMFT calculations presented above, the residual resistivity at T = due to random 
scattering from the incoherent i sites is too large, ~ ImVLcm (Fig.3). If the i states 
become coherent, the b electron scattering and the consequent resistivity would then 
vanish at T = 0, and is nonzero only if static disorder is present. This can lead to a 
metallic state with a small residual resistivity or to an Anderson localized insulating 
state depending on the amount of disorder. 

We believe that the giant isotope effect observed in manganites [HI] is another 
dramatic consequence of the £ — b hybridization scale t± and its exponential depen- 
dence on the inverse square root of the isotopic mass. For it can add, via (con- 
ventional) double exchange, to the ferromagnetic T c an amount roughly given by 
AT C ~ afii(D* /ks) — (87a) K for ng = 0.7 where a is less than one. Because of the 
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Huang Rhys factor in D*, which depends exponentially on square root of the isotopic 
mass, we find that AT c (0 16 ) ~ (1.33) AT c (0 18 ) so that for a ~ (1/2) the difference 
T c (0 16 ) - T c (0 18 ) ~ 15 K, close to the observed value(H]. 

In addition to the neglect of intersite coherence effects discussed above, we have 
also neglected in our theory spatial correlations between the local electron densities 
nn and ribi, as well as between the local angles Q{. This may be adequate for the 
phases which correspond to homogenous orbital liquids. However, many phenomena 
in manganites such as short range order, long range charge/orbital order (or orbital 
solidification) , various types of antiferromagnetism |3J E] and mesoscale structures 
jH] depend on spatial correlations. To treat these we need to add to eq. (1) a number 
of more complicated longer range coulomb, anharmonic, steric, elastic, magneto- 
elastic, etc. interactions that couple spin, orbital and lattice degrees of freedom 
to each other, and to strain, ion size mismatch, disorder, etc. as appropriate. A 
self consistent determination of long or short range order in 5i = {fin — fibi) (a 
new 'internal' variable) [48J , 6i and Si, and their effects on the b electron dynamics 
for different x and T, in the presence of such interactions can lead to a complete 
description of manganites including the above phenomena, and others such as first 
order transitions, two phase co-existence, etc.. We hope to discuss these elsewhere. 

A significant general question raised by our work is that of the detailed nature of 
adiabatic to non adiabatic crossover as g increases, and the conditions for the coex- 
istence of adiabatic and antiadiabatic states with exponentially separated dynamical 
timescales. We have shown that the latter leads to new phenomena, argued that 
this happens in manganites because of the large g, and have developed the conse- 
quences in a simple model assuming this separation. The results closely correspond 
with a wide variety of observations. It would be of great interest to explore such a 
crossover and time-scale or energy-scale separation experimentally and theoretically 
in the many systems such as organic solids, transition metal oxides and chemical 
(molecular) systems, which have degenerate orbitals and strong symmetry breaking 
JT couplings. 
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